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Abstract: The use of model membranes is currently part of the daily workflow for many biochemical and biophysical disciplines. 
These membranes are used to analyze the behavior of small substances, to simulate transport processes, to study the structure of 
macromolecules or for illustrative purposes. But, how can these membrane structures be generated? 

This mini review discusses a number of ways to obtain these structures. First, the problem will be formulated as the Membrane 
Packing Problem. It will be shown that the theoretical problem of placing proteins and lipids onto a membrane area differ 
significantly. Thus, two sub-problems will be defined and discussed. Then, different — partly historical — membrane modeling 
methods will be introduced. And finally, membrane modeling tools will be evaluated which are able to semi-automatically generate 
these model membranes and thus, drastically accelerate and simplify the membrane generation process. The mini review concludes 
with advice about which tool is appropriate for which application case. 

Mini Review Article 



Introduction 

The basis for many biophysical and biochemical in silico 
experiments is a model membrane. Experiments may be used to 
assemble protein-membrane systems to verify and simulate the 
structure of proteins or to theoretically evaluate values derived by wet 
lab experiments (e.g., the area per lipid, the membrane thickness, the 
viscosity with respect to temperature changes and the transition 
between different phases). The behavior of small organic and 
anorganic substances can be evaluated as well as the permeability, 
active and passive transport processes. And of course these model 
membranes can also be used for illustrative purposes. 

These structures have to be generated using fragmentary 
knowledge about the composition of biological membranes since it is 
not possible to acquire the complete membrane structure with 
presently used microscopy or spectroscopy methods. As there is no 
standard procedure, there are many different ways to acquire these 
initial structures. Basically, three different methods can be 
distinguished to generate a model membrane. 

(1) Structure Recycling: First, they may be ^'recycled''. This means, 
they can be taken from ones own previous in silico experiments or 
downloaded from websites and other repositories. They can also be 
taken from the supplement of published material. The advantage is 
that the membrane layer is already equilibrated. However, it is not 
always possible to find a simulated membrane showing the 
composition needed for a computational experiment. 

(2) Manual Structure Generation: These structures can be generated 
by scripting or by composing them using a text editor. For example. 
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there are some software tools available which, although they were not 
originally developed to generate membranes, may be used in 
conjunction with a script to generate membranes. Related approaches 
will be summarized under the chapter Membrane Modeling Methods. 

(3) (Semi-) Automated Structure Generation: In the past few years 
new software tools were developed which were optimized for the 
modeling of membranes..These tools will be introduced towards the 
end of this mini-review. 

This mini review explicitly cannot provide an extensive summary 
of all tools and methods available today. But to the knowledge of the 
author, this publication covers the most important ones. This 
discussion will primarily focus on the creation of rectangular 
membrane patches. This work concludes with suggestions for the 
appropriate tool for a given application case. 

But first, the theoretical background of membrane modeling has 
to be discussed. This background was previously introduced in our 
publication from 201 1 [I]. 

From Nature to Model 

Nature is our role-model. Since the Fluid Mosaic Model was 
introduced in the 1970s, many different observations have approved 
and extended this theoretical membrane model [2]. It is known that 
the membrane is not a rigid structure. It is a highly dynamic 
culmination of different specific biological entities such as proteins 
and lipids. These entities are in continuous fast-motion. Lipids change 
their positions and structure. Some of them are changing from the 
external side of a membrane bilayer to the internal side [3]. Proteins 
bind with lipids or other proteins located inside the membrane. 
Substances are transported from one side of the membrane to the 
other. And of course, while the elements of the membrane are 
changing, the shape of the membrane changes too. Lipid rafts float 
like small islands through the membrane sea [4] and membrane caves 
are formed, ready to enclose cargo, separate from the membrane and 
transport the cargo through the cytosol. 
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Figure 1. The result of the MOE Membrane Grid script as shown in Jmol. 



Very complex problems arise when the membrane modeler takes 
all these different aspects into account. How is it possible to simulate 
all the biological, chemical and physical aspects of all those membrane 
components? The physicist might wish to use quantum mechanical 
approaches, but then it would be nearly impossible to simulate the 
structure of a single protein. Therefore, approaches exist combining 
molecular dynamic and quantum mechanical approaches [5]. 
Although pure quantum mechanical approaches are applied to a 
problem, it is still not possible to prove that the resulting in silico 
structures absolutely represent their real state in vivo, 

A number of different approaches exist to simulate membrane 
structures, trying to reproduce experimentally derived values in the wet 
lab. The initial guess could be that these simulation approaches could 



be used to generate an initial membrane model. This approach will be 
shortly discussed as the self-assembly method. The problem is that the 
intended structure will often not be achieved. If, for example, the 
membrane modeler wants to generate a rectangular membrane patch, a 
system could be created containing water and a number of lipids 
randomly placed in a cubic box. If a Molecular Dynamics (MD) 
simulation is now applied to this system, the lipids will usually tend to 
form micelles instead of generating rectangular structures. Here is a 
major dilemma because the rectangular membrane patch is a fragment 
which usually does not occur in nature. But it is applied to many 
different application cases in order to examine a small fragment of 
reality as the molecular simulation of a complete cell membrane lies 
beyond the scope of present computer technology. 
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An alternate approach is to generate a stiff starting structure using 
different techniques. Usually, in the next step, it will need to be proven 
that the resulting in silico structure is stable over time by applying a 
molecular simulation. In this case, the experimentalist will need force 
fields as well as the coordinate information of the lipids. The force 
fields describe the physico chemical behavior of lipids, including 
torsion angles and atomic forces. This information is only available for 
a small amount of lipids. However, all actual simulation environments 
such as CHARMM, GROMACS, or NAMD, provide force fields 
compatible with lipid simulations [6—8]. 

Packing Problems 

If the process of creating a membrane is broken down into its 
most basic task(s), it can be understood as a packing problem. These 
are known from a wide range of economical and ecological 
applications. One example is the logistic sector, where the task is to 
optimize the relation between space and cargo [9]. Packing problems 
are known to be NP-hard [10], and solving them is a complex task. 

Following the definition of our previous publication, the issue is 
described as the Membrane Packing Problem (MPP) [ //. Two sub- 
problems have to be distinguished. First, the Lipid Packing Problem 
(LPP) and second the Protein Packing Problem (PPP). 

• (LPP): Lipid Packing Problem: How can one or multiple lipid 
models be optimally packed onto a monolayer or bilayer? 

• (PPP): Protein Packing Problem: How can one or multiple 
protein models be optimally packed onto a monolayer or 
bilayer? 

• (MPP): Membrane Packing Problem: How can one or 
multiple lipid and protein models be optimally packed 
together onto a monolayer or bilayer with respect to LPP and 
PPP? 

In the following chapters it will be shown that it is important to 
make the distinction as the packing of lipids onto a membrane area 
differs significantly from the insertion or attachment of a membrane 
protein. 

Membrane Modeling Methods 

A number of traditional membrane modeling methods should be 
discussed: the Grid-based, the Replacement, the Insertion and the 
Self-assembly method. As previously mentioned, these — partly very 
time-consuming — methods were especially relevant in the past. Then 
there were no software tools freely available which were optimized for 
the membrane generation process. But there are a number of relevant 
PPP methods listed in the following sections which are not integrated 
into semi-automatic Membrane Modeling Tools discussed towards the 
end of this manuscript. 

Often, molecular modeling packages had to be used which 
provided a scripting interface. The user had to write a script which 
created an initial membrane model, as discussed in the next section. 
An alternative way was to write small customized tools. Often they 
were never published by the authors and are therefore not freely 
available. All methods discussed here are partly used to create pure 
membranes as well as protein/membrane complexes. 

Grid-based Methods 

The idea of the Grid-based method is quite simple: the molecules 
are placed along a virtual grid, representing the horizontal layer 
located between the two layers of a bilayer membrane. If 64 lipid 



models should be placed at one layer, a virtual grid with a resolution 
of 8x8 is generated, and a molecule is placed at each intersecting point 
of this grid. 

A typical script should be shortly explained. It was re- 
implemented by the author of this mini review based on the method 
published by Kriiger and Fischer [16]. For this purpose, MOE 
(Molecular Operating Environment) was used. A number of 
alternative tools could also be applied. MOE is a commercial and 
well-established modeling tool which is used, for example, in drug 
discovery [17]. MOE is not originally intended to create membrane 
structures. But it is possible to do so by writing a script, using the 
scripting vector language (SVL). The applied script implements two 
loops. The first loop operates along the Z-axis, the second loop along 
the X-axis. A defined number of lipids are placed along each axis with 
respect to the user-defined distance. The user may also choose if a 
monolayer or bilayer should be created. The resulting DPPC lipid 
with ten lipids per row and a distance of 10 A can be seen in Figure I. 

By examining the result in Figure I it is obvious that it is not able 
to meet biological structural expectations, because the grid is too 
regular. The resulting structure is very vague and ignores geometrical 
or chemophysical aspects. All molecules in the resulting structure 
feature the same distance from the neighboring lipids. Therefore, in 
the best case, a number of minimization steps wiU be needed to 
generate a membrane model compatible with molecular simulations. 
But in the worst scenario, the gaps between the lipids might be too 
large. This state will usually result in a membrane, where the water 
molecules, added during the solvation process, fill the bilayers' gaps. 
During the simulation, the water-filled gaps will expand and in this 
way they will usually destroy the bilayer structure. Of course, the 
structure of the membrane in Figure I could be improved by choosing 
initial grid values meeting experimental expectations concerning the 
area per lipid. It will be densely packed but it can be surmised that 
the atoms of the lipids will collide or even intersect. An add-on 
minimization might rectify the situation, although some simulation 
packages may have extreme difficulty with these intersecting 
structures. However, to prevent as many collisions as possible, it is a 
good idea to use lipid models with stretched-out tails. It can be 
summarized that this is a quick method to generate an initial structure 
for those researchers familiar with a molecular modeling environment 
providing scripting interfaces. But in the next step it might be 
troublesome because the energy minimization might take too long or 
could even fail. In both cases it would be a better approach to use a 
packing which tries to avoid collisions. 

In addition, this membrane is only homogeneous. The MOE 
Membrane Grid script could be extended or the replacement method 
discussed in the following section could be applied if different types 
of lipids have to be inserted. 

The membrane modeling tools discussed in the next chapter could 
be practical for the researcher who has not previously written a script 
or who is unfamiliar with a similar tool. 

An interesting method combining a Grid-based method with 
protein placement was discussed by Kandt, Ash and Tieleman. They 
applied a widely-spaced grid of lipids, similar to the one shown in 
Figure I, while placing a protein in the grid's center. Then, they 
shrank the grid around the protein, moving the lipids closer to the 
edge of the protein, until the intended area per lipid was reached 
[12]. 

Replacement Methods 

Initially, Pastor et al. described the Replacement method [18], 
Different publications on this process are summed up in the following 
section: 



Volume No: 5, Issue: 6, February 201 3, e20l3020l4 



Computational and Structural Biotechnology Journal | www.csbj.org 



Membrane Modeling 



If the final membrane structure should contain a single protein or 
peptide, this molecule may be used as the starting structure 
[19,20]. During the following steps (Step 2 to 7), this molecule 
defines the center of the membrane and it remains fixed at its 
initial position. 

The starting point of the lipid membrane is a single lipid model 
which is used during the generation process. But instead of using 
the regular atomic structure of this lipid, it is represented by a 
sphere. This object may be the representative for a headgroup 
[18], a Lennard-Jones sphere of 4.8 A radius corresponding to its 
average cross-sectional area [19], a van der Waals radius sphere 
[20] or spherical beads with an approximated size of the polar 
headgroups [21]. 

Now, the spheres representing the lipids are distributed onto a 
plane (by omitting protein in the center from step 0, if present). 
The placement of the lipids may follow a regular hexagonal [18] 
or octagonal arrangement [19], or they are just randomly 
distributed onto the plane [20,21]. 

This step might not be necessary if homogeneous bilayers are 
constructed [18]. But the packing of spheres representing the 
lipids might have to be optimized if heterogeneous bilayers are 
constructed. For this purpose, different approaches may be used. 
This might be an MD simulation [19,20] as well as a shrinking 
process using the steepest descent minimization algorithm 
[21,22]. During the different repositioning processes, the vertical 
positions of the sphere usually remain fixed. But it may also be 
that the vertical movement is possible within a short range [20]. 

Until now the spheres were utilized to substitute the atomic 
structure of the lipids. Now that the distribution process is 
finished, these spheres have to be substituted by the original lipid 
structures. These are randomly chosen from a trajectory, which is a 
library containing a limited number of different conformations of 
a lipid chain [18]. These trajectories have been previously created 
by simulations, often taken from already published and verified 
data [23,24]. If the orientation of the lipids is not uniform, the 
lipids may also be rotated randomly in a range of 0 to 360° in the 
X/Z-plane and tilted in a range from 0 to 45° with respect to the 
bilayer which is normal during the substitution process [18]. 
In this context, the placing process does not usually take collision 
detection into account. For this reason, there will be steric 
collisions inside the generated membrane model if it is imported 
into a molecular simulation environment. This may be problematic 
if a molecular simulation should be invoked afterwards. 

For this purpose, different methods are now applied which are all 
similar or equal to minimization processes. During the energy 
minimization of the membrane structure, the radii of the atoms 
grow from I to 100% [18]. Another approach is that the lipid 
structures (and optionally also the protein or peptide in the 
center) remain stiff while they are systematically horizontally 
moved and tilted [19]. Sometimes this method can also be 
combined with an increasing process of the atom radii [20]. 
Sometimes the Steepest Descent (SD) minimization algorithm can 
also be used for this task [21]. 

These optimization methods reduce the bad contacts but usually 
they do not disappear (for example, a decrease of factor 2 is 
observed [19]). To solve the dilemma of the remaining collisions, 
a number of atomistic energy minimizations has to be applied. 

Now, the model membrane is ready for the equilibration process. 
Therefore, all constraints, like the fixation to a position in the 



X/Z plane, are removed and the molecules are able to move and 
interact freely. 

8. Finally, the initial structure is complete and the molecular 
simulation, e.g., an MD simulation, can be started. 

Insertion Methods 

The Insertion method is often used to integrate larger molecules 
such as proteins or peptides into a previously generated bilayer (PPP). 
Therefore, the subsequent — partly historical — examples roughly 
discuss protein insertion methods with focus on the surrounding 
bilayer. In contrast with the previous method (Replacement method. 
Step 0), first the membrane is generated, and then the protein is 
inserted into the surrounding layer. A short overview of methods will 
follow. 

1. A combination of the insertion method with the replacement 
method was discussed by Shen, Bassolino and Stouch. The base 
was a previously-published starting structure of a well- equilibrated 
bilayer membrane. At each layer, two lipid models were removed 
from the center. But the hole inside the bilayer was not large 
enough to accommodate the peptide. To increase the radius, weak 
cylindrical repulsive forces were appended to the hole during the 
subsequent minimization and MD simulation process. As soon as 
the radius of the whole was large enough, the peptide was 
inserted. But because the atoms of the lipid and those of the 
peptide were still colliding at some positions, the system was 
minimized again. Finally, after all collisions were eliminated, the 
equilibration and MD simulation could follow [25]. 

2. A classical approach was discussed by Kandt, Ash and Tieleman. 
The protein is inserted into a pre-equilibrated bilayer and 
overlapping lipids are removed based on cut-off lipid-protein 
distances. The aforementioned publication lists a number of 
approaches with this strategy. As the authors state, the main 
disadvantage of this approach is that the distances between the 
remaining lipids and the protein are irregular. Therefore, 
difficulties might arise when starting the simulation because the 
initial simulation box size might already shrink during the first 
simulation steps [12]. 

3. A slightly different approach was applied by Tieleman and 
Berendsen. Initially, a bilayer containing 64 lipids on both sides 
was created using a grid-like positioning and random rotating 
method. After the short solvation and simulation process (using 
periodic boundary conditions) the bilayer was multiplied by 4, 
resulting in a membrane containing 256 lipids at each side. Then, 
two methods were used to insert a protein into the membrane. 

a. First, the protein was inserted into the membrane without 
removing intersecting lipids after which an abstract two- 
dimensional grid was laid on top of the membrane area. 
Each quadratic area was checked for the coexistence of the 
protein with lipids. All lipids affecting the quadratic area 
were removed. The advantage of this approach is that the 
resulting membrane did not exhibit any atomic collisions. 
On the downside, it showed large distances between the 
lipids and protein. 

b. Alternatively, Tieleman and Berendsen used the previously- 
discussed approach but with the exception of allowing 
minor intersections. After the initial structure was created, 
a number of energy minimization attempts combined with 
force-related investigations and the removal of some 
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intersecting lipids were needed. It was more complicated 
to achieve then the first approach but the resulting 
structure did not show large gaps between lipids and 
proteins [26]. 

4. Faraldo- Gomez, Smith and Sansom began with the generation of 
a solvated and equilibrated lipid bilayer. Next, the solvent- 
accessible surface area of the protein was computed. The resulting 
shape was used in the following process as the protein's 
placeholder. Then, all lipids were removed if their headgroup P 
atom was intersecting a cylinder- shaped region where the protein 
should be placed. In the next step the protein was inserted into 
the bilayer and only some lipid tails extended into its region. The 
small number of intersections were eliminated during the 
following MD simulations [27]. 

5. Instead of placing a protein, it is also possible to place other 
specific molecules into a virtual membrane area. One example was 
discussed by Deleu et al. In this case, a surfactin was fixed in the 
membrane's center. Then, a new lipid entered the membrane area 
and moved along the X-axis in small increments less than an 
Angstrom. At each step, the lipid slightly rotated, changed its 
vertical position and tilted. Meanwhile, the energy of each 
composition was calculated and the one with the lowest energy 
was kept using a tool called HYPERMATRIX. Then, the next 
lipid was added and went through the aforementioned procedure 
[28,29]. 

6. Yesylevskyy introduced a plugin called ProtSqueeze. It was 
initially developed to be used in conjunction with VMD and the 
source code could be used with virtually any other program 
supporting plugins or scripts. The workflow of ProtSqueeze will 
be shortly summarized [30]. 

a. A pre-equilibrated bilayer has to be used and a protein has 
to be pre-aligned into the bilayer. For this purpose 
ProtSqueeze is able to use OPM structures for the pre- 
alignment [15]. Alternatively, ProtSqueeze provides a 
number of scripts, which try to place the protein in 
relation to the bilayer based on its residues defining the 
two membrane-water interfaces. Afterwards, all lipids are 
deleted which intersect the protein based on the steric 
clashes. The user may define a percentage value which 
varies the grade of intersection. 

b. Now, the squeezing process starts. At this stage, the 
structure of the protein shrinks until all steric clashes are 
eliminated. A special method is used to take any possible 
tilts of the protein into account. 

c. Finally, the shrunken structure of the protein has to be 
expanded again. At each sub-step, the protein is slightly 
stretched, followed by an energy minimization of the 
surrounding lipids to remove possible new overlaps. Step 3 
is repeated until the protein's original size is reached, of 
course avoiding atomic collisions. 

A very similar approach is today part of the GROMACS 
simulation suite. It is called g_membed [13]. 

7. Jo, Kim and Im developed an insertion method in conjunction 
with a web interface (to be introduced later in the Membrane 
Modeling Tool section) [31]. A hole is created inside a pre- 
equilibrated bilayer by applying repulsive radial forces around the 
center of the protein. For this purpose, lipid bilayers were 
generated, equilibrated, and collected in a library which contains 



membranes of various sizes and hole radii. Based on the size of 
the protein, an appropriate bilayer is chosen and the position of 
the protein-surrounding lipids is adjusted based on the repulsive 
radial forces. 

8. Staritzbichler et al. introduced a methodology called GRIFFIN, 
which improves the protein insertion process. Lipids and water 
molecules of a pre-equilibrated system are removed corresponding 
to the volume of the protein. The protein is placed in the space 
created. However, still non-protein atoms intersect the protein. 
Then, a grid-based implicit force field is used during a subsequent 
Molecular Dynamics simulation. Atoms of non-protein molecules 
lying inside the protein's volume are confronted with an outward 
force which moves them outside the protein's volume, whereas the 
external non-protein atoms are stirred by van der Waals and 
electrostatic forces. The big advantage is that the overall volume 
and density of the simulation box is preserved and the 
surrounding non-protein molecules are smoothed to the surface 
of the protein [14]. 

9. Another option placing a protein inside a membrane is to acquire 
protein pre-alignment information from specific databases. 
PDB_TM and OPM provide information which enables the 
alignment of proteins contained in the PDB database with respect 
to the center of a bilayer. For this purpose, PDB_TM offers a 
transformation matrix which has to be applied to the structure of 
the original PDB model. OPM provides a different approach; it 
attaches dummy atom layers to the PDB model, indicating both 
sides of the membrane's bilayer. Both approaches use special 
algorithms to compute the protein's membrane alignment. 
Moreover, OPM offers a service which tries to attach the dummy 
layers to custom PDB files uploaded to the OPM server. Of 
course, these methods can be combined with the aforementioned 
Insertion methods. Both pre-alignment approaches can also be 
semi-automatically used in conjunction with the 
CELLmicrocosmos MembraneEditor discussed in the Membrane 
Modeling Tools section [15,32]. 

Self-assembly Methods 

Another method is the self-assembly of membrane layers. It is 
theoretically a simple experiment for users of molecular simulation 
programs. A number of lipids are randomly placed in a virtual box 
using simple manual methods such as a text editor or molecular 
modeling program. These lipids are then immersed in water and then 
a molecular simulation is started. After a relatively short period 
(usually less than a nanosecond), the self- assembling process starts. 
The lipids often will try to form a vesicle. The head groups are 
hydrophilic and will orient themselves towards the water interface. As 
the tails are hydrophobic, they are directed towards the inner center of 
the vesicle. 

Of course, this is only a very rough description and the success of 
this experiment depends on several factors, the applied force field, the 
duration, the temperature and pressure settings, etc. But there are three 
major problems. (I) The user is not able to exactly define the shape of 
the resulting structure. So it will be quite complicated to generate a 
bilayer model. (2) The challenging requirement on computational 
resources. If a well-established MD simulation package such as 
GROMACS is used, it is usually a good idea to use a computer or 
GPU cluster [7]. (3) Preparing such an in silico experiment is very 
time consuming if the user is not familiar with molecular simulation 
approaches. 

Of course, there is also a big advantage: the resulting structure is 
in a perfect state for further simulations. 
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Geometry-based Methods 

To prevent the lipid's distribution resulting in an overly artificial 
structure, it might be useful to begin by applying geometry-based 
methods during the generation process. This means that atom-based 
or shape-based collision detection is taken into account. The 
Geometry-based method is presented here as a serious alternative to 
the simple but fast grid-based method and the more realistic but 
complex methods described previously. 

In our previous publication, we defined the assembly of a 
membrane layer as a Two-and-a-half-dimensional Knapsack Problem 
(2.5D-KP) [I]. The basics have to be explained to describe this term. 
For this purpose, the Lipid Packing Problem has to be examined from 
a computational perspective. One well-known computer science- 
related packing problem is the Knapsack Problem (KP), which is 
formally defined as follows [33]: 
n 

maximize ^ ^/ , 
/= 1 

n 

subjectto constrain ^W-X- < X • G {0,1,...,C.} (1) 

i=l 

There are n kinds of items x with values of Vi and a weight of Wu 
And the overall weight of the given container is restricted to Here, 
the term container defines the distribution area. The number of copies 
Xi is restricted by the maximum value c which is a characteristic 
criterion for the Bounded KP (BKP), Therefore, the knapsack's area 
(for a 2D-KP) or volume (for a 3D-KP) is limited, in contrast with 
the Unbounded KP (UKP), where a is defined as infinity. The values 
V and w are positive integer values /' in the context of the LPP 

The original fluid mosaic membrane model already proposed the 
membrane as a two-dimensional liquid arrangement [2]. Now the 
question is, if this definition is sufficient in the context of LPP? 

The knapsack's weight Wi is often used to describe dimensionality. 
An extension of the KP is the so-called Multidimensional KP where 
knapsacks may feature multiple dimensions [33]. In the context of 
describing rectangular bilayers, the two-dimensional [9,34] or 
geometric\iS\ KPcIslss is appropriate (2D-KP). 

The basic KP is restricted to one container. But there is also 
another extension of the KP called the Multiple KP (MKP) which 
includes multiple related containers [36]. In the case of the LPP, a 
monolayer can be defined as one container and a bilayer as two 
containers. In contrast to the MKP, these two containers of a bilayer 
are not related. For example, in the case of the MembraneEditor, the 
modeling process strictly focuses on the geometric properties of the 
molecules, refusing chemophysical interactions as well as periodic 
boundary conditions. Therefore, both layers are regarded as being 
strictly divided by the hydrophobic core in the center of the bilayer. 
Of course it is known that in vivo interactions between both layers 
steadily occur, for example in case of cholesterol flip-flops [24]. 

Now, the different variables of Equation I have to be assigned to 
the LPP. The container represents one membrane layer whereas the 
different items are the molecules. For the regular LPP, many items of 
relatively few differing figures or shapes exist [36]. The items are the 
lipid models and their frequency is defined by the lipid ratio. In the 
case of the LPP, the weight Wi is represented by the width (along the 
X-axis) and length (along the z-axis) of the molecule. This might be 
confusing for the experimentalist, because a biological lipid has 
molecular weight. But this weight is not relevant for a strictly 
geometry-based approach. But of course, the height and width are 
important, because a larger lipid needs more space in the knapsack. 
The value Vi is one, because each molecule is counted as one and 



added to the overall number of lipids. If different lipid types are part 
of the layer (for example, 80% phospholipid and 20% cholesterol), Ci 
is restricted by the lipid percentages and BKP would be the 
appropriate class. In this case, the lipid will not be placed in a free 
space if the actual percent value for the according lipid type is equal 
or larger than the lipid percentages. By contrast, UKP is the 
appropriate definition if the according layer contains only one lipid 
type. Then, Ci is infinity. 

Many classical KP solutions are not applicable to LPP, because the 
item's orientation is also fixed or it is only altered in 90° steps [9,34]. 
This is not possible if the items are represented by lipids. Packing 
problems such as KP are NP-hard [37,38]. Various alternative non- 
exhaustive solutions exist, for example: approximate algorithms, greedy 
algorithms, and heuristic algorithms [33]. A single good solution is 
sufficient for many application cases in contrast to an exhaustive one 
which is usually not computable during an appropriate period of time. 

Finally, the term 2D-KP should be reevaluated. For strictly two- 
dimensional problems — for example, the distribution of boxes or 
spheres on a two-dimensional area — this definition is correct. Even if 
the shape of the lipid on a two-dimensional plane is irregular — which 
is almost always the case — it is still a two-dimensional packing 
problem. But in the case of distributing three-dimensional objects like 
molecules, the collision detection should also operate in three 
dimensions, whereas the movement of the lipid is restricted to the 
two-dimensional plane. 

Therefore, the term 2.5D-KP was introduced in our previous 
work [I]. It was shown, that the 2.5D-KP is the appropriate problem 
definition for the geometrical packing of lipids on a plane. 

Discussion 

The last two sections showed that the borders between the 
different methods are floating. They are often combined. For example: 
(I) The lipid's layer is created using a Grid-based method and then a 
replacement or insertion method is applied. (2) Lipids are initially 
substituted by simple structures but the proteins remain as atomic 
structures, or, inversely (3) the proteins are temporarily represented by 
shapes and the lipids are handled as atomic structures. 

The traditional generation processes are naturally quite complex. 
They are often not straight- forward. The previous examples showed 
that often multiple minimization attempts are needed to generate a 
membrane model usable for the following simulations steps or during 
the minimization of the molecules' radii (which have to increase to 
perform the alignment and so forth). 

Another problematic aspect of most traditional Membrane 
Modeling Methods is that the generation protocols found in many 
articles are vague and therefore the results are not reproducible. Of 
course, it is an ongoing debate in science if all methods described in a 
publication must be reproducible. The following example is one 
among a thousand which illustrate this aspect. (In no way should the 
quality of the following publication be doubted as this sentence is 
only shown for illustrative purposes.) 

''The configurations were assembled as a set of rigid units, with each 
GA or DMPC ( with their primary waters ) being translated and 
rotated in a systematic search for an optimum packing'' [39]. 

But how do the traditional Membrane Modeling Methods apply 
to the KP? The initial problem is the same; an area has to be filled 
with lipids. Moreover, the replacement method utilizing spheres as 
lipid-placeholders can be interpreted as a typical 2D -UKP. The 
knapsack is unbounded if the spheres are all of the same type. It is a 
two-dimensional packing problem, because the movement of the 
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spheres is restricted to the X/Z-plane and the packing can be reduced 
to the two-dimensional area where the radii of the spheres are at a 
maximum size. After replacing the spheres with the lipid's atomic 
structure, different methods are used to remove collisions. Methods 
like the minimization and equilibration, which are normally used for 
the replacement methods, are not typical packing problems, since 
chemophysical methods are applied here which are not restricted to 
geometrical issues. Still, the Replacement methods could be seen as a 
packing problem but with different new criteria. But these criteria are 
not restricted to the shape. They apply to every single atom during the 
minimization processes. Every atom can move in different directions. 
But regular packing problems, generally and more particularly the KP, 
take stable shapes or outer boundaries into account. The LPP is 
defined for stiff lipid structures and the KP optimization takes 
advantage of the fact that the lipid's structure does not change during 
the distribution process. The criterion **many items of relatively few 
different figures or shapes'' would not be met in this case [36]. 

Membrane Modeling Tools 

Five tools will be briefly introduced providing capabilities to 
generate membranes without the need for computing complex 
molecular interactions on a local computer. In addition, these tools do 
not assume chemophysical expert knowledge and the regular computer 
user is able to adapt the tool's workflow in an appropriate time frame. 
Each of the tools discussed here must be able to export structures in a 
format which is compatible with MD packages. And the regular 
export format which is supported by all introduced tools is the 
Protein Data Bank (PDB) format [40]. Although the name of this 
format suggests that is only appropriate for proteins, nowadays it is 
the standard format for all applications working with three- 
dimensional molecular structures. It is, for example, possible to obtain 
lipid structures based on molecular assemblies from the PDB 
repository using the websites HIC-UP or Ligand Expo [41—43]: 

http://xray.bmc.uu.se/hicup/ 
http:/ /ligand-expo.rcsb.org 

It has to be emphasized that these are tools solving the Membrane 
Packing Problem. In the section discussing the Membrane Modeling 
Methods, there were also tools listed which explicitly solve the Protein 
Packing Problem (ProtSqueeze, g_membed„ and GRIFFIN, for 
example). 

ChemSW® ChemSite Pro ® 

Chemsite Pro® (Version 7) from the company ChemSW® is the 
second commercial tool besides MOE. This package is also to model 
and simulate molecular structures in addition to other functionalities 
[44]. It provides an intuitive approach for modeling membranes. This 
program also comes packaged with the MolSuite™ which provides 
quantum mechanical simulations and a database containing different 
molecular structures. 

The module Build Lipid Membrane of ChemSite Pro is 
interesting in the context of this review. It provides a number of 
simple options to create a lipid membrane. First, the user can decide if 
a symmetric or asymmetric monolayer or bilayer should be created. 
Then, the width of the quadratic -shaped membrane is defined by the 
number of lipids along each edge multiplied by the user-defined 
distance between each lipid. It is not possible to generate membranes 
with different side lengths. Obviously, this is a Grid-based approach. 

It is possible to apply a simple solvation to the membrane by 
defining the number of water layers at each side of the bilayer as the 
addition of water layers is often needed for MD simulations. 



ChemSite Pro places parallel rows of water layers beneath the heads 
of the highest lipids. Therefore, distances will appear between small 
lipids like cholesterol and the water layer. The resulting membrane 
models will be problematic for a sophisticated simulation, because the 
layer might collapse when the water molecules penetrate the bilayer 
and they are found everywhere in the system, even between the two 
lipid layers in the hydrophobic area. 

Now, the lipids have to be defined. Theoretically an unlimited 
number of lipid types can be added, each represented by a name and a 
percent value. In case the membrane is not symmetrical and a bilayer 
was chosen, the values for each lipid type have to be defined for both 
sides and it is possible to generate asymmetric bilayers. Finally, the 
Build'' button has to be pressed and the membrane is created. If the 
contained lipids are charged, counter ions are automatically added to 
maintain a state of equilibrium. 

The disadvantages of the grid-like structure have been previously 
discussed. Another problem is that the resulting membrane is not a 
mixed bilayer because the different lipid types are placed iteratively; 
the starting rows are filled with lipid type I until the percentages are 
reached, and then the following rows are filled with lipid type 2 and so 
forth. 

ChemSite uses a native library format for lipids, called lib. It 
contains the coordinate information and some further information 
regarding charges. It is not directly possible to use PDB files to 
generate these membranes. It is, however, a simple task to translate 
PDB to lib files by using ChemSite Pro, to re-import them and to 
integrate them into the membrane. The lipid library which comes with 
ChemSite Pro contains six different types, four phospholipids, one 
phosphatic acid and one cholesterol. 

Moreover, ChemSite Pro includes a nice tool to create a large 
number of custom lipids, the Lipid Builder 

But there are also some disadvantages. The program exhibits low 
performance. For example, if a membrane with 100 lipids is explored, 
juddering occurs unless a powerful computer is used. Moreover, the 
layout of the GUI is outdated and there seems to be no community 
supporting or using this tool. Because of the very simple generation 
method, the resulting membranes exhibit an extremely high energy 
state. Therefore, a very long minimization run is needed before the 
MD production run can be started. This is problematic as it it not 
possible to use a cluster system for the minimization process with 
ChemSite Pro and it must be computed on a desktop system. 

VMD Memhrane-Plugin 

YMD (Visual Molecular Dynamics) is a popular program to 
visualize, analyze, and edit molecular structures and simulation files 
[45]. It is also possible to simulate membranes by connecting external 
packages like NAMD [6]. VMD also supports plugins, such as the 
''Membrane" plug-in which was developed by Ilya Balabin, on-line 
since 2006 [46]. It is the simplest tool discussed here, but the fact 
that it is still part of the actual VMD release 1. 9. 1 shows that it is 
accessible to a large community. The simple syntax of the membrane 
plugin's command line tool is as follows: 

membrane -1 <lipid_name> -x <size_in_X> -y <size_in_Y> 
{-o <output_prefix>} 

The lipid type is defined behind the -1, where there are only two 
options: POPC (a phosphatidylcholine) and POPE 
(phosphatidylethanolamine). Only bilayers containing one single lipid 
type are possible. The lipids feature outstretched tails to improve the 
packing quality. The size of the rectangular membrane is defined by 
the variables -x and -y. The variable -o defines the prefix for the two 
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membrane output files. These are a PDB file (containing the 
coordinate information) and a PSF file (containing force-field-specific 
structural information of the membrane model). 

The membrane is generated by combining hexagonal patches of 
shortly equilibrated solvated POPC or POPE bilayers. The solvation 
is based on the Solvate plugin from Grubmiiller and Groll which 
generates an irregularly-shaped disordered water layer [47]. This 
resulting structure can be directly used via VMD for MD simulations. 

The tool, packaged with VMD, is found at: 
http://www.ks.uiuc.edu/Research/ vmd/ 

CHARMM'GUI Membrane Builder 

CHARMM ( Chemistry at HARvard Macromolecular Mechanics ) 
is a highly established program package providing different tools and 
force fields for the minimization and simulation of different 
molecular structures [8,23]. Unlike the free MD tools such as 
GROMACS, the user has to obtain academic or commercial licenses. 

The CHARMM-GUI initiative is quite new and it may be used 
free of charge. Here, the version 1.2 is discussed. As denoted by the 
name, the target group is mainly users of the CHARMM package. 
CHARMM-GUI provides a set of tools which should simplify the 
process of generating starting structures for MD simulations [48]. 
The usage is quite simple because CHARMM-GUI can be directly 
accessed via the website [49]: http:/ / www.charmm-gui.org 

For this mini review, only the Membrane Builder will be discussed. 
It is part of a collection of tools found in the website section called 
'Input Generator''. At the moment, it is only possible to generate 
bilayers; monolayers are not possible. For homogeneous bilayers, six 
lipid types are available (DMPC, DPPC, DOPC, POPC, DLPE and 
POPE) and for heterogeneous bilayers, 32 lipid types can be used 
(incl. cholesterol and different sphingolipids). Each lipid type in this 
library is represented by 2,000 different conformations extracted from 
MD simulated bilayers [50]. It is possible to use a quadratic as well as 
hexagonal shape. 

Now, the generation process for ''Membrane Only Systems*' will 
be discussed. It is subdivided into five different steps. 

Step 0) In the preliminary step, the conformation of the bilayer is 
defined, including the shape, the thickness of the water layer, the 
initial estimate for the rectangular size and whether or not ratios or 
numbers should be used as a reference point. All available lipids — 
depending on the type of membrane — are listed and the user can 
choose the absolute or the percentaged distribution values for each 
bilayer side. It is also possible to edit the initial guess for the surface 
area for each lipid type. There is a button which can be used to 
evaluate the predicted properties of the membrane model which will 
be generated as follows: 

Step I) The initial files are generated and the properties of the 
files are now shown. For the assembling process the user can choose 
between the Insertion method (only available in homogeneous 
bilayers, see the corresponding Membrane Modeling Method of Jo, 
Kim and Im) or the Replacement method, from Woolf and Roux. 
This latter method places pseudo atoms onto the membrane area 
[19]. After the initial distribution has finished, the pseudo atoms are 
replaced by molecular structures of the aforementioned lipid library. 

The insertion method depends on a library of pre- equilibrated 
bilayers which are restricted to a maximum size of 90x90 A^, 
containing maximal 256 lipids. For the Membrane Builder, the rapid 
Insertion method should be used for proteins of regular and 
cylindrical shape and for all other cases the Replacement method. If 
charged molecules were found, the program is able to inhibit counter 
ions and the user may choose between different ion types [31]. 

Step 2) shows the previously- created membrane components; the 



ions, a newly generated water box and the lipid-holding layer. It is now 
possible to download and/or to look at a preview of the PDB 
structure of the lipid layer using an integrated Java applet called 
MarvinSpace (Version 5.4) from ChemAxon Ltd. Note; this tool 
seems to have problems with larger membrane files. 

Step 3) provides information and editable options regarding the 
equilibration process, the temperature, the ensemble type and the 
surface tension. 

Step 4) enables the user to download the resulting files which can 
be used for further equilibration and simulation with CHARMM. 
These files are also compatible with NAMD simulations which can be 
used with VMD (discussed in the previous section). 

The workflow and the results provided by CHARMM-GUI are 
quite powerful, particularly for users of the commercial CHARMM 
software package. It also provides output files compatible with 
NAMD, so users of this software package might be interested in using 
this tool. Although initially planned, CHARMM-GUI does not 
currently provide simulation files compatible with GROMACS [23]. 

An important and unique feature of the Membrane Builder in 
comparison to the other tools discussed here is the computation of 
lipid's surface areas by using Voronoi tessellation based on Pandit et 
al. [51]. These estimations are appropriate criteria for judging the 
space occupied by a lipid model. In the case of heterogeneous 
membrane layers, this method is more accurate than just computing 
the average area per lipid, because the shape among different lipid 
types strongly differs. 

Although the generation process is well-described in publications, 
the source code is not available. Users of Open Source software, such 
as GROMACS, prefer to use software providing access to all data of 
the simulation process. Therefore, it might be undesireable to use a 
module in the modeling and simulation process which shows 
similarities to a black box. 

In terms of reproducibility there are also some drawbacks. First, 
the program is only available via Internet and reproduction of the 
results is only possible as long as the website is accessible and the 
provided methods are not changed by the developers. While the use of 
the Membrane Builder is presently free of charge, further simulation 
using the CHARMM simulation package requires a paid license. 

Because it is not possible to import custom lipid models, many 
MD packages like GROMACS will probably require the additional 
manipulation of the PDB membrane created by the Membrane 
Builder by using the Replacement method for example. 

The generation methods do not take collision-detection of 
covalent or van der Waals radii into account. Therefore, collisions 
occurring after applying the Replacement method have to be 
eliminated by more expensive minimization and other special methods 
described in the original publication of the Membrane Builder [50]. 

Packmol 

In the collection of tools now presented, Packmol is the only tool 
which is exclusively available as a command line tool. Therefore, its 
handling will not be easily accessible to a number of users. But it has 
to be discussed because it is quite powerful as it is able to generate 
molecular structures of different shapes such as spheres, cylinders, 
planes and/or boxes by applying geometry-based methods. As with 
most other tools discussed here, it supports homogeneous as well as 
heterogeneous membrane layers. The old version of Packmol has 
already been used for many publications since 2003 [52]. Published in 
2009, the new version supports a more effective distribution 
algorithm combined with the ability to parallelize the computation 
[53]. 
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It uses geometry-based methods implemented using a special 
packing algorithm called GENCAN [54]. It is able to solve the 
previously introduced 2.5D-KP, but it is also able to generate 
membrane structures with complex three-dimensional shapes like 
vesicles. This tool is particularly powerful if small molecules need to 
be packed into three-dimensional volumes defined of various shape. 

However, Packmol is also often used for generating rectangular 
homogeneous membranes. Moreover, there are also a small number of 
publications found which generate heterogeneous membranes with 
Packmol, e.g.. Hall et al. [55]. For this purpose, the user exactly 
calculates the absolute number of lipids which will be placed inside 
the membrane layer. The lipid models are manipulated to be 
compatible with Packmol. Since the position and orientation along 
the vertical axis of the different lipid models should match they are 
pre-aligned. 

Packmol is invoked by a customized script. This script contains all 
information about how the membrane has to be constructed. For new 
users it might be a good idea to search for an existing script and 
customize it. The positioning of lipids can be defined via constraints. 
It is even possible to restrict the vertical position of the headgroups to 
be placed beneath a predefined plane. This feature is quite powerful 
but requires the user to have an exact idea of the coordinate space. 
Therefore, mathematical abstract thinking is a prerequisite for the 
correct use of Packmol. 

The packing process integrated in Packmol may result in a 
number of problems during the generation of a membrane. If the 
structures are too complex, the packing optimum might never be 
reached. To find out which layer sizes and lipid ratios are compatible, 
the user often has to find out by trial and error. In addition, the 
resulting structure has to be equilibrated if MD packages like 
GROMACS, CHARMM or NAMD are used. Packmol does not 
offer any visualization capabilities. But of course, for this purpose 
external tools like Jmol can be used [56]. The program is located at: 
http://www.ime.unicamp.br/~martinez/packmol/ 

CELLmicro cosmos 2.2 MembraneEditor ( CmME) 

The CELLmicrocosmos 2.2 MembmneEdior ( CmME) is the 
newest of the tools represented here. It was developed for the fast and 
user-friendly modeling of membranes. It is the only tool strictly 
applying geometry-based methods to the LPP In contrast to 
CHARMM-GUI Membrane Builder, it does not support the 
equilibration of the membrane. The philosophy of the CmME is that 
the equilibration should be done directly with the tool of (^oic^, like 
GROMACS or NAMD. In our publication from 201 1 it was shown 
that these methods produce good results by using GROMACS [I]. 
Moreover, we are currently developing the GMX-Plugin which should 
be a bridge between CmME and GROMACS. 

A number of Lipid Packing Algorithms (LPA) provided by 
CmME solves the 2.5D-KP. Rectangular, homogeneous or 
heterogeneous monolayers as well as bilayers can be computed using 
different plug-in algorithms. It is possible to use percentages to define 
the composition of the membrane as well as absolute values. Micro 
domains/lipid rafts can be drawn by hand or shaped as rectangles or 
ellipses. It is also possible to generate stacked membranes using 
CmME. 

Every valid PDB file can be imported and defined as a lipids or 
protein. Lipids are automatically aligned perpendicular to the lipid 
layer. The semi-automatic placement of proteins which solve the PPP 
by using PDB_TM and OPM was already mentioned in the 
Membrane Modeling Methods section as well as in the original 
publication of the MembraneEditor [I]. There are two ways of 



combining a protein and a membrane. First, it is possible to use the 
Insertion method. In this case, the lipids are distributed on the plane, 
the protein is inserted and then the intersecting lipids are removed. 
Problems for transmembrane proteins arising from this procedure 
were already discussed in the Insertion method section. In the other 
case, first the protein/s is/are placed onto the empty membrane area 
and then the lipids are distributed around the protein using the Lipid 
Packing Algorithms. For transmembrane proteins, this is the 
procedure of c^oict. For extrinsic proteins, the Insertion procedure 
should be used. 

CmME features WYSIWYG (What You See Is What You Get). 
The lipids are visually represented by shapes and their distribution is 
also directly visualized. Therefore, associations with the replacement 
method are possible. But this does not completely apply to the 
approaches of CmME, because the collision detection may take the 
visual shapes and/or the atomic structures into account. Moreover, 
the structure of the atomic structure can also be directly visualized 
with CmME. But it is advised to use Jmol — which is bundled with 
CmME — for this purpose. 

To judge the packing density, CmME computes the average overall 
surface area. In contrast to CHARMM-GUI Membrane Builder, this 
method is not based on Voronoi tessellation and is therefore less 
accurate. But it provides initial insight into the quality of the LPAs. 
Moreover, we are currently developing a program which provides 
extensive analysis functionalities for lipid membranes by using Voronoi 
tessellation. 

It was previously mentioned that Packmol supports the generation 
of vesicles. For CmME, there is also an initial version of a Vesicle 
Builder available at the same website. Here CmME can be directly 
downloaded and started as Java Web Start application: 
http:/ / Cm2.CELLmicrocosmos.org/ 

Conclusions 

Five different Membrane Modeling Method categories were 
introduced: the Grid-based methods, the Replacement methods, the 
Insertion methods, the Self-assembly methods and the Geometry- 
based methods. By examining these methods it was shown that the 
distinction of Membrane Packing Problems into I) Protein packing 
Problems and 2) Lipid Packing Problems is reasonable, because in 
many cases only one of these problems was addressed. Furthermore, if 
proteins and lipids are combined in a system, the packing processes of 
both molecule types usually differ significantly. 

Then, a number of Membrane Modeling Tools were discussed 
which are all able to handle the Membrane Packing Problem and 
which implement at least one of the aforementioned Membrane 
Modeling Methods. 

MOE® was used to present a Grid-based method which can also 
be applied by using the scripting capability of other tools. This 
approach is interesting for all those users who, I) already know their 
program of c^oict quite well and, 2) are familiar with the scripting 
language. For all other users it will be more efficient to use one of the 
following approaches. 

ChemSite Pro® from ChemSW® might be good choice for those 
people using the MolSuite™ developed by the same company. 
However, this grid-based approach is quite simple and there are gaps 
between the placed lipids. Moreover, different lipid types in 
heterogeneous membranes are not mixed. Therefore, model 
membranes which should be simulated — especially those with external 
packages — should be modeled using one of the following tools. 
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Membrane Modeling Tools Overview 

Table 1. The complete comparison of all membrane modeling tools discussed in this work (Charmm GUI MB: Charmm GUI 
Membrane Builder; VMD MP: VMD Membrane Plugin; GMX: The Gromacs-Plugin for the MembraneEditor; VB: The Vesicle Builder 
Plugin for the MembraneEditor). 



Category 


Feature 


Charmm 
GUI MB 


Chem 
Site Pro 


CmMRi 


MOE+ 
Script 


Packmol 


VMD MP 




Standalone 




X 


X 


X 




X 




Command line tool 










X 






Web service 














Availability 


Website 


X 












Web Start 






X 










Source Code 






X 




X 


X 




Licenses 

(educational license fee) 




commer-cial 


GPL3 


commer-cial 


GPL3 


UIUC Open 
Source 


Computational Acceleration 


Multi-threading 














Parallelization 










X 




Direct PDB Database 




X 




X 


X 




















GRO 












X 


Formats 


native format 




X 


X 


X 








PDB 


X 


X 


X 


X 


X 


X 




Grid-based 




X 




X 




X 




Replacement 


X 




X 








Membrane Modeling Methods 


Insertion 


X 




X 










Self-assembly 
















Geometry-based 






X 




X 






Free 






X 




X 






Hexagonal Shape 


X 












Membrane Shapes 


Quadratic (X=Z) 


X 


X 


X 


X 


X 






Rectangular (X?^Z) 






X 


X 


X 


X 




Vesicles 






(VB) 




X 






Lipid library (>2 lipid types) 


X 


X 


X 








Libraries 
















Lipid library compatible to 
MD 


X 


X 


(GMX) 










Average Surface Area per Lipid 






X 








Lipid Area Computation 
















Surface Area by Voronoi 
Tessellation 


X 














Asymmetry 


X 


X 


X 




X 






Atom-based Molecule Editor 




X 


X 


X 




X 




Automatic Alignment of 






X 










imported Lipids 
















Bilayers 


X 


X 


X 


X 


X 


X 




Collision-Detection 






X 




X 






Counter ions support 


A 












Features 


Lipid Packing Density 


X 




X 










Heterogeneity 


X 


(X) 


X 




X 






Lipid Ratio: Absolute 


X 




X 


X 


X 


X 




Lipid Ratio: Relative 


X 


X 


X 










Mono-layers 




X 


X 


X 


X 






Multi-layers 






X 




X 






Protocol 






X 










Raft Support 






X 




X 







Reproducibility 






X 










Water Layer Build 


X 


X 


(GMX) 


X 


X 






WYSIWYG 




X 


X 


X 




X 




Linux 


X 




X 


X 


X 


X 


operating System 


Mac OS X 


X 




X 


X 


X 


X 




Wmdows XP, VIS 1 A, 7 


X 


X 


X 


X 


X 


X 


Pipelines 


Ext. Simulation Package 


X 


X 


(GMX) 


X 




X 


Ext. Visualization Package 


X 


X 


(Jmol) 






X 


Programming 


Scripting 






X 


X 


X 


X 


Scripting at Runtime 






X 






X 




Equilibration 


(X) 




(GMX) 


X 




X 


Simulation 


Minimization 


(X) 


X 


(GMX) 


X 


X 


X 




Simulation 




X 


(GMX) 


X 




X 




Graph Visualizations 






X 


X 








Live Distribution Visualization 






X 








Visualization 


Atomic Structure View 


X 


X 


X 


X 




X 


Secondary Structure View 






Qmol) 


X 




X 




Ray-tracing 




X 












Runtime Graphs 






X 










Stereo Support 






X 


X 




X 



^GMX = GMX-Plugin (in development); Jmol = Jmol library (directly included in CmME); VB = VesicleBuilder-Plugin (in development) 



The YMD Membmne-Plugin is the right c^oict for modeling a 
homogeneous POPC or POPE bilayer membrane. The biggest 
advantage is that it is part of the popular VMD and nearly every 
member of the modeling community will have this tool on the 
computer. This plugin is the first c^oict if the membrane should be 
simulated with NAMD and there is no need for other lipid types. 

CHARMM-GUI Membrane Builder is a very powerful website- 
based tool. It applies the Replacement as well as the Insertion method 
(only for homogeneous bilayers) to solve the Lipid Packing Problem. 
It usually will be the first dvoict for all users of the CHARMM 
simulation package. It provides a large number of regularly -used lipid 
types. However, it is not possible to import custom lipid types. 
Another small drawback is that the shapes of the membranes have to 
be quadratic. For users of other simulation packages — especially 
Open Source tools — one should think twice about the use of the 
Membrane Builder because the source code is not available and the 
tool is optimized for use with CHARMM. But the additional option 
to generate NAMD -compatible files is a step in the right direction. 
The Membrane Builder is also well-described in its original 
publication [50]. 

Finally, the two Lipid Packing approaches contributed by the 
Open Source community have to be summarized. 

Packmol is a very powerful and versatile package which is able to 
generate molecular structures of many different shapes using 
geometry-based methods. One application case is the generation of 
lipid layers. In contrast with the following, it is a command-line tool 
providing no visualization and is invoked with a script. The writing of 
this script needs the user to have a clear picture of the geometrical 
properties of the molecules used and the resulting structure. The 
lipids have to be pre-aligned and the composition of the membrane 
has to be defined by constraints. Packmol directly tries to minimize 
the resulting structure by using its custom methods. In any case, it will 
be the first c^oic^ if the user needs to generate structures consisting 
of very small molecules such as water, which should be three- 
dimensionally distributed [53]. 



CELLmicrocosmos 2.2 MembraneEditor ( CmME) was developed 
from the start to generate heterogeneous membrane patches in a fast 
and easy way by using geometry-based methods. For a fast modeling 
and evaluation of a membrane composition it should be the first 
<do\z^. Of course, if chemophysical properties of a membrane should 
be evaluated, the resulting structures should be minimized, 
equilibrated and simulated by external, well-established simulation 
environments like GROMACS. We recently showed that structures 
generated with CmME are compatible with MD simulations of 
GROMACS. A first version of a Vesicle Builder is also available. The 
Java Web Start version of this tool enables quick installation of this 
tool [I]. 

As the title of this mini review indicates, the detailed explanation 
and evaluation of the Membrane Packing Problem lies beyond the 
scope of this exposition and the complex methodologies are only 
roughly described. Moreover, this short review does not claim to list 
all known Membrane Packing approaches. 
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